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In this paper we consider a stochastic process that may experience random reset events which 
bring suddenly the system to the starting value and analyze the relevant statistical magnitudes. 
We focus our attention on monotonous continuous-time random walks with a constant drift: the 
process increases between the reset events, either by the effect of the random jumps, or by the action 
of the deterministic drift. As a result of all these combined factors interesting properties emerge, 
like the existence — for any drift strength — of a stationary transition probability density function, or 
the faculty of the model to reproduce power-law-like behavior. General formulas for two extreme 
statistics, the survival probability and the mean exit time, are also derived. To corroborate in 
an independent way the results of the paper, Monte Carlo methods were used. These numerical 
estimations are in full agreement with the analytical predictions. 
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I. INTRODUCTION 

Stochastic processes that incorporate some kind of 
mechanism which may instantaneously bring the mathe- 
matical system to an absorbing state are not new P, : 
The motivation may arise from many real-life systems 
which are threatened by the possibility of a sudden and 
severe decimation, usually with an external origin, e.g., 
the extinction of biological species due to some catas- 
trophe, bankruptcy of industrial and financial firms or a 
major electric-power outage. 

Those situations where such a catastrophic event re- 
sults in the disappearance of the system have attracted 
considerable interest and are well studied in the litera- 
ture, both from a mathematical and physical perspec- 
tive. Mathematical models that encompass this possi- 
bility arc fairly developed: For diffusion processes (like 
the Brownian motion) the possibility of extinction has 
been well described and several killing mechanisms have 
been considered in the literature. The introduction of 
an absorbing boundary at the origin, say, is the simplest 
of them. In this situation the system vanishes once it 
hits the given level. However, the latter does not reflect 
all possible physical situations. Neither it does exhaust 
the possible behaviors at the boundary, which with more 
generality may display a combination of absorption, re- 
flection, elastic, sticky and jump behavior — see for 
the classification of boundaries for diffusion processes. A 
different mechanism, where boundary behavior plays no 
role, assumes that in any time interval (t, t + At] the 
system has a probability of extinction V{x)At where x 
is the actual position and V{x) the killing rate. Unlike 
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what happened previously, in this setting extinction oc- 
curs in a sudden way, with no forewarning and may take 
place at any location x and time i 0, Q . 

In this paper we are interested in a related but never- 
theless fairly different situation where the system is not 
extinguished but rather is "reborn from ashes." In such 
a case the (exogenous) event represents a reset in the 
system activity: the reboot of a computer, the neural 
spiking activity, or the "back-to-square-one" order of cer- 
tain board games. The hallmark of such a system is the 
possibility of a random restart, in opposition to random 
disappearance. Hence this novel behavior cannot be cap- 
tured by considering alternative sorts of boundaries, like 
reflecting or sticky ones. Actually, here we are interested 
in the opposite case where the triggering of resets is to- 
tally independent of the underlying system. 

In spite of its obvious interest and significance, it ap- 
pears that this situation has, until very recently, been 
virtually ignored in the literature. As far as we know, 
this possibility was first considered in the 1999 seminal 
work by Manrubia and Zanctte Q • This paper analyzed 
a discrete-time stochastic multiplicative process with a 
finite event space ^ namely a Markov chain endowed with 
a restart mechanism. The second exception is due to 
Evans and Majumdar 0,0] who in 2011, apparently un- 
aware of the previous work of Manrubia and Zanette, 
considered the diffusive evolution of a particle that may 
be randomly reset to its initial position. A remarkable 
conclusion, found both in Q and Q , was that the shut-off 
mechanism induces an algebraic decay for the tails of the 
corresponding probability density functions. Power laws, 
i.e., "heavy-tailed" distributions, have been repeatedly 
observed in natural plwsical phenomena and advocated 
by several authors [^-|lOl| but, despite this fact, power- 
law tails are by no means ubiquitous: most mathematical 
models ignore them in favor of Gaussian decay — see [ll| 
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for a critical analysis of this dichotomy. 

Following these ideas, in this paper we analyze the in- 
troduction of a restart mechanism in the framework of 
continuous-time random walks (CTRWs) with a constant 
drift and study the consequences that this carries. Thus, 
the dynamics of a such system involves three elements, 
namely the random jumps, the drift term and finally the 
resets. 

If we deprive the system of the CTRW part but con- 
serve the drift and the random reset mechanism we re- 
cover the classical shot-noise process, a well known model 
to describe sudden light emission with current induc- 
tion 

Conversely, pure CTRWs correspond to having only 
the first of these elements. The interest of these pro- 
cesses in Physics has been advocated specially by Mon- 
troU and Weiss [3 ■ Statistical Physics driftlcss 
CTRWs have been used to describe transport in disor- 
dered media fl6l-l20l , earthquake modeling [21I . [22| , ran- 
dom networks |23| . self-organized criticality in granular 
systems 24 1, electron tunneling [2^, and distribution of 
matter in the universe [26j . just to name a few. But 
CTRWs have been also used in the modeling of social sys- 
tems, e.g ., by gi ving a tick-by-tick description of financial 
markets [27l - l30| . A comprehensive review of CTRW ap- 
plications in Finance and Economics can be found in l3ll . 
For a more general perspective on the issue see |32h34| |. 

While not so popular in Physics, the CTRW with a 
drift incorporated plays a fundamental role to model the 
cash flow at an insurance company [35l - [38j . Recently, we 
have shown that energy dissipation of optical beams in a 
sclf-dcfocusing medium like an optical fiber with random 
inhomogeneities is also described in terms of a CTRW 
with a drift [H, H^l ; see also [4l| . Typically, in such a sit- 
uation signal's losses may be significant and may require 
the incorporation of periodic all-optical amplifiers [42| 
to upgrade the signal strength to the initial amplitude. 
Thus, in this context, the exogenous reset mechanism 
proposed here might be incorporated to the model in a 
natural way. 

The paper is organized as follows: In Sec. |TT] we in- 
troduce the process under study, a CTRW with drift 
on which a restart mechanism has been super-imposed. 
We suppose that resets are set off with no forewarning 
and independently of the value of the underlying CTRW. 
The resulting process is both rich and realistic enough to 
deserve important attention. In Sec. IIIII we obtain an 
explicit analytic expression for the the transition prob- 
ability or propagator of the system. It turns out that 
the full propagator can be related in a relatively simple 
way with the rcsct-frcc one. We also find that as long 
as the reset mechanism is present, the asymptotic be- 
havior (in cither time or space) shows some appealing 
features, like, in particular, the existence of a stationary 
distribution. We pinpoint conditions that guarantee the 
existence of power-law tails and characterize precisely the 
relevant exponent. Section IIVI is devoted to the analy- 
sis of the properties of two extreme-value statistics, the 



survival probability and the mean exit time. We give a 
closed formula for these statistics in terms of a Laplace 
transform. Conclusions are drawn in Sec. |Vl where fu- 
ture perspectives are also sketched. We have left for the 
appendices some technical mathematical derivations. 



II. THE PROCESS 

Our starting point is a random process X{t) whose 
dynamics consists of a superposition of a constant rate 
motion along with sudden random increments J„, the 
jumps, at random instants of time i„. More concretely, 
if To and Ti are given times we suppose that for To < 
t <Ti, X(t) can be expressed as 

00 

X(t)=rx (t-To)+ Jn0{t-tn)9{t^-To), (1) 

n— — cxD 

where F > is the constant velocity, and 9{-) is the 
Heaviside step function, that takes values 9{u) = 1 for 
u > and zero otherwise. 

Furthermore. t„ and J„ are the jump times and jump 
sizes while Tn = tn — tn-i are the time intervals between 
two consecutive jumps, the "waiting times." We assume 
that J„ and r„ define sequences of independent and iden- 
tically distributed (i.i.d.) random variables, described by 
the probability density functions (PDFs) h{-) and tp{-), 
respectively: 

h{u)du = P {u < J,i < u + du} , 
ip{T)dT = ¥{t <T„ <T -hdr} . 

As we are interested in the analysis of monotonous pro- 
cesses, we will further assume that h{u) = for u < 
(whereas iP{t) = {01 t < by definition.) 

Up to the time Ti Eq. ^ defines simply a CTRW 
with a constant drift F like that considered in [4l|. The 
difference with previous works arises at i = Ti , when the 
process experiences a reset event,^ i.e., we declare that 
X{Ti)^0 and that for Ti < t < T2, X{t) is given by 

00 

X{t)=Tx{t-Ti)+ Jn9{t-tn)9{tn-Ti). (2) 

n— — oo 

The same definition applies to X{t) at t ~ T2,T^,..., 
i.e., XiTm) = 0; where Tm is an increasing sequence of 
random times. Thus in general for Tm < i < Tm+i, 
m = 0, 1, 2, . . ., we have 

00 

X{t) = Tx{t-T„,)+ Y Jn9{t-t,,)9{tn-T,n). (3) 



^ Our mathematical model assumes that both jumps and reset 
events do suddenly modify the state of the process. While in 
some cases this instantaneous-response hypothesis will be an ac- 
curate description of the system, in others could only be consid- 
ered as an acceptable proxy of reality. 
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Note how past Ti the initial evolution of X{t), ef. 
Eq. ([T|) , is modified in an exogenous way. The novelty of 
this paper is the introduction of this mechanism within 
the CTRW framework, which brings forward the possi- 
bility that the system be shut down and then restarted, 
see Fig. [T] 
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Figure 1. (Color online) Sample path of the processes X{t) 
and Y{t). The solid (red) line represents a possible realization 
of the monotonous process X(t) which grows linearly between 
the jump times t„ and returns to the origin at the reset times 
T„. The dashed (blue) line is a geometrical version of X(t), 

y(t) = e-^w. 

As regards the sequence of random resets we as- 
sume that the inter-event intervals 5^ = Tm — Tm-i are 
a set of i.i.d. random variables characterized by the PDF 

(j){T)dT EE P {r < ^„ < r -f dr} . 

We also assume that the set {S^m} is independent of {t„} 
and { Jn} and hence that resets occur independently of the 
value of the underlying CTRW . 

The process X{t) thus modified is intrinsically different 
to the original one, and, generically, is not a Markov 
process. In fact, it is not either a compound renewal 
process: On the one hand, if i„ is a given jump time, 
the evolution of X{t) after that time depends on X(i„) 
but also on the previous reset time T*{tn) = Tk, where 
Tk < tn < ?fe_|_i, for some k. On the other hand, if T„i 
is a given reset time, the evolution of X{t) after that 
time depends as well on the last jump time, which we 
denote by t*{T„i) = ti, with ti < Tm < U+i for some 
value of I. In other words, in spite of the formal resem- 
blance of Eqs. dl]) and (0), X{Tq + r) and X{Ti + r) have 



In spite of this fact, one could argue that the long-time limit of 
the new system might follow by using an infinite combination of 
systems without resets. While in the analysis of some properties 
of the process this can be true, see Eq. I I14I I below, it is unclear 
if this approach could be substantiated to produce concrete ana- 
lytic results in the general case. We are grateful to an anonymous 
referee for this observation. 



different laws and therefore one cannot resort to the use 
of renewal equations in the analysis of this process for 
a general choice of ■(/'(■) and as is commonly done 
when studying CTRWs. 

In this paper we take the natural and usual choice that 
the waiting-time and reset densities ip{-) and are 
exponential functions 

^(r) - Ae-^^ (4) 
0(r) = Ae-A^ (5) 

where and A~^ are parameters of the distributions 
that respectively denote the mean value of t„ and 3^ra- 
Note that with this choice X{t) is time- homogeneous and 
Markovian. Indeed, if no resets are present X{t) re- 
duces to the well known compound Poisson process with 
drift which satisfies the above properties. (The inter- 
ested reader can find the proof of this fact in App. |B1 
along with several global properties of this process.) 

It is interesting to point out that in many cases 
the physically measurable process Y{t) is not the given 
by Eq. ^ but is rather constructed by multiplica- 
tive — instead of additive — increments. To accommodate 
this interpretation we simply consider that X{t) is the 
(aggregated) rate at which the physical observable^ Y{i)^ 
varies. In such a case X{t) and Y{t) arc related either 

by 

y(0=loe^(*), (6) 

or by 

r(0 -^oe-^^*). (7) 

This is the interpretation that we take in the rest of this 
paper where X{t), defined by Eq. ([3]), is mostly consid- 
ered for computational convenience, whereas the magni- 
tude of interest is Y{t). Let us enumerate a few examples: 

In Finance, Eq. ^ may represent the evolution of a 
savings account with a minimum balance level Kq (lot us 
say, one cent) and a continuously compounded interest 
rate F, wherein sudden injections of capital are made by 
the owners at random times i„ , but which does not allow 
for partial withdrawals. 

In Biology, Eq. ^ may model the changes in the popu- 
lation of a species in some geographically-limited region, 
assuming that it has an effective growth rate F and is 
subjected to recurrent immigration and punctual mas- 
sive emigrations: e.g., the case of many birds, swallows 
or flamingos among others. In this case the reset event 
docs not imply the extinction of the species as a whole, 
but the return to the autochthonous level of population 

In an optical context, the evolution of the energy E{t) 
of a signal in an optical fiber with random impurities is 
also described by a related model, as we we have recently 
shown [3§,|43l. Here E[t) = Eq e-^(*), where X{t) is de- 
scribed by Eq. ([T]), F accounts for the constant damping 
rate of the fiber, while jump points f„ signal the presence 
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of random impurities in the fiber due to manufacture er- 
rors. As a result of these combined effects, the initial 
signal will be degraded exponentially in time. The intro- 
duction of an ulterior reset mechanism, as we do here, 
could be used to describe the possibility that random 
(digital) amplifiers — at which the signal is amplified to 
the initial value, i.e., E(Tn) = Eq — are incorporated to 
the fiber to make up for the losses due to the damping 
and impurities. 

Finally, we note that when A = 0, namely when only 
resets but no jumps are allowed, we have that Eq. (O 
leads to 

oo 

Y{t) = YoX ^ e-^^'-^-'>9{t-T^)e{T^+i-t). 

m— — oo 

Therefore, Y(t) is the piece-wise exponential function 
Y{t) = Foe-^^*-^") for < t < T,„+i. This is the 
shot-noise process where shots are triggered with Pois- 
sonian rate and suffer an exponential decay after hit- 
ting the detector. Shot noise has been used widely in 
Physics to describe sudden electron emission with 

induced currents, and Cherenkov radiation. 

Thus, the system considered here is a natural gener- 
alization of classical physical models. In the rest of the 
paper we study the most relevant statistical magnitudes. 



Equation ^ intends to emphasize that, in principle, 
the propagator depends on both the starting point xq 
and the initial time to- However, the process must be 
time-homogeneous (since both the jump and the reset 
mechanisms are Poissonian), and therefore 

p{x, t; xo, to) = p{x, t ~ to; xq, 0) = p{x, r; xq). 

By contrast, one does not expect the propagator to be 
translationally invariant because resets take the process 
to a fixed point, thereby breaking the spatial symme- 
try. This implies that the dependence on xq cannot be 
removed. 

Let us derive first the integral equation that governs 
the evolution of p(a;, r; xq). Recall that and t (with 
to < t) are, respectively, the actual and a future given 
time instant, and that t = t — to is the associated time 
interval. In the following we will denote by t' and r" the 
time interval up to the first reset event and, respectively, 
the first jump — note that r' and r" are random variables 
while T is just a number. The equation for the transition 
PDF can be built up by considering the three possible 
and mutually exclusive scenarios that appear depending 
on the relative values of the time interval up to the first 
reset event, r', the time interval up to the first jump, r", 
and r: 



III. THE TRANSITION PDF 

We start with the analysis of the propagator, the tran- 
sition probability density function of the process, 

Px{x, t;xo,to)dx = V {x < X{t) < x + dx\ X{to) = xq} , 

(8) 

where t and to are arbitrary |j with t > to- 

Note how in the definition of pxix,t] xo,to), the "X" 
subscript indicates the random variable to which the 
transition PDF is associated; similarly, we denote by 
Pviy, t] yo, to) the corresponding probability for the phys- 
ical process Y{t) = e^^^\ However, since py follows 
trivially from px , 

PY{y,t;yo,to) = -pxilogy,t;\ogyo;to), (9) 

y 

we will mainly focus our attention on px{x,t; xo,to) in 
the sequel and, to simplify the notation, drop the "X" 
subscript when there is no room for confusion. 



(i) There is neither a reset nor a jump in the time in- 
terval T, i.e., r' > r and r" > r. In this case 
X{t)^xo + Tt. 

(ii) There is at least one reset in the time interval, and 
the first one takes place before any jump has oc- 
curred, t' < T and t" > t' . In this case the transi- 
tion PDF after the reset must be p{x, r — r'; 0) since 
system must reach x, starting from zero, in a lapse 
of time T ~ t' . 

(iii) There is at least one jump in the time interval (to, t], 
namely t" < r. This first jump (of size u) takes 
place before the first reset event, r' > t". Right 
after the jump we have X{t") = xo + Ft" + u, and 
then the propagator is p{x, t ~ t"; xq + Ft" + u). 

In view of all this p{x,T;xo) must satisfy the following 
renewal equation: 



We note that in the previous section we used to to denote one of 
the jump times. This abuse of notation is not merely incidental. 
When the process is not time-homogeneous the renewal equations 



have a simpler structure if the time origin coincides with a jump. 
* We set lo = 1 hereafter for computational convenience. In prac- 
tical terms, this can always be achieved by a suitable choice of 
the physical units. 
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dr"Ae 



dT"\e-^^5{x-XQ-VT) 



dT"Xe-^^ p(.T,T-r';0) 



-At 



duh[u)p[x, T 



+ / duh{u) / dr"Ae-(^+^)^"p(x,T-T";a;o + rr" 



';a;o + rr" + u) 
t';0) 



(10) 



The standard procedure for solving integral equations 
like (Hn]) is to resort to the use of some integral transfor- 
mation, either the Laplace transform, the Fourier trans- 
form, or a combination of them. Here, since x and r are 
positive variables, the natural choice is to consider the 
Laplace transform in both arguments: 



p(w,s;a;o) 



dre 



dxp{x, r; Xo)e 



where from now on the hat will denote the Laplace trans- 
form with respect to the time variable, and the tilde will 
denote the Laplace transform with respect to the posi- 
tion variable. Typically upon transformation the inte- 
gral equation yields a purely algebraic equation whose 
solution can be (straightforwardly) found. In the present 
case, however, the double Laplace transform simplifies 
but does not solve the posed problem. 



p(uj,s-xo) 



1 



\ + A + Tuj + s 

A + A + s 

dr"Ae-(^+^+-^)^" 



duh{u)p{uj,s;xQ + TT" + u), (11) 



since p{uj^ s; xq) satisfies a new (although simpler) inte- 
gral equation. 

Further progress in the resolution of Eq. (|lip demands 
a careful analysis of the dependence of p{x,T]Xo) on xq 
or, in other words, determining which is the interplay be- 
tween xo, X and r. In App. |2] we prove that the solution 
to this integral equation is 



p{uj,s;xo) 



where 



1 



1 - h{Lo) 



A + Fw- 



- s 
(12) 



duh{u)e 



The propagator p{x,t]Xq) then follows by inversion 
of the Laplace transform of Eq. (|12|) . Letting A — ?> 
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one recovers pq(x,t]Xq), the propagator of a reset- free, 
monotonous CTRW with a drift, the process that follows 
if we let Ti ^ c» in Eq. ([T]), 



Po(w, s; Xo) = lim p(uj, s; xq) = — 



1 - h(uj) 



■ Tuj - 



(13) 

It is noteworthy that the opposite construction can also 
be carried out, namely p(a;, r; xq) can be written in terms 
oi Pf){x,T]XQ). In physical space this relationship reads 
(see App. El 



p{x,t;xo) 



[ dT*Ae-^^'po{x,T*;0) 
Jo 



+ po{x,T;xo)e 



-At 



(14) 



Equation p4p is more informative from a structural 
point of view than Eq. (fT2|) since it relates both propaga- 
tors. Actually the proof given in App. [B]on the Marko- 
vianity of X{t) relies heavily on this fact. Nevertheless 
we find it more convenient to use Eq. (jl2p as starting 
point for obtaining p{x, r; xq). In particular, one can go 
one step further and take the inverse Laplace transform 
in the s variable of p(a;, s; xq) to obtain p(uj, r; xq): 



5{uj,t;xo) 



A 


1 


- e 


(A[l- 




A+rw)r 




A 


1 - 


h{uj) 


+ A 


+ Fw 




-(A 


l-/i(w)]+A+ 


raj)r 



(15) 



An ulterior inverse Laplace transform in the uj variable 
could be given once the functional form of /i(cj) is chosen. 
Before examining further the general properties of the 
propagator we first consider a concrete example. 



A. Exponential jumps 

Let us exemplify the general results in Eqs. (|12p - P^ 
with the following choice for the jump density 



h{u) =7e-T", {u > 0), 



or equivalently 



h{uj) 



7 



7 + a; 



(16) 



(17) 
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where 7 > is a parameter. This election is both plausi- 
ble and suitable, but even in this case the calculations are 
substantial. Since the purpose of this section is mainly 
illustrative we first consider the simpler case F = 0. 
Under these assumptions Eq. reads 



We close this section with some remarks concerning 
the stationary density for the general case F > 0. The 
inclusion of a non-vanishing drift constant blurs some- 
what the neat result in Eq. (|22p with the emergence of 
two different exponents 



A 

-+e- 
s 



(A A s) w + (A + s) 7 ' 
(18) 

from which one can readily compute the inverse 
Laplace transform: 



p{x,s;xo) = 



1 



A + A + s 
A 7A 



— S{x) + S{x — xq) 
s 

(A+3)7X 

_e A+A+s 



s {X + A + sY 

2^ p-XTiTT7(— -0)/ 



(A + A + sY 



{x - Xo). 



(19) 



To completely recover the density one must, in addition, 
to perform the inverse Laplace transform in the s vari- 
able, which is in turn more intricate. After some calcu- 
lations (see App. [C]) we finally obtain 



p{x,t;xo) 



A 



-Six) 



-(A+A)t 



X + A 

S{x ~ xo)e-^^+^^^ 



A 



2y/x¥ja 



g-(A+A)T-7(2:~a;o) 



7Ar 

x — Xq 



X h (2y^XT-/{x - Xo)) 0{x ~ Xo), 



(20) 



where /i(-) is the modified Bessel function of order 1. 

Formula ([20)1 implies the existence of a stationary 
PDF, which follows taking the limit r — >■ 00 and com- 
puting the integral that appears in the third term. Al- 
ternatively, one finds in a straightforward way from (jl9p 
that 



p(x) = lim p(x,t;xq) = lim sp(x, s; xq) 
- — s->-0 



r— >-oo 

A 
A + A 



Six) 



7A _A2SL 
-g A + A 



A-f A 



(21) 



Thus, the stationary distribution of the physically 
relevant quantity Yir) = e^^''^ follows from Eqs. © 
and (I2T]) 



pviy) 



1 A 



5i\ogy) 



yX + A 



7A A-y log y 



A-l- A 
7AA 



1 



(A + A)^y 



2 ■ 



with 



7A 
A + A' 



(22) 



(23) 



Pviy) 



A 



F(q;+ — «_) 



a+ — 7 7 ■ 



il+a- 



(24) 



with 



A + A + F7 , 

a+ = ± 

2F 



y^(A + A + F7)^ 



4F7A 



2F 



(25) 



Since one has a+ > 7 > a_ > 0, the bigger parameter 
a+ dominates the behavior of py (y) for y ~ 1 whereas 
the smaller one, a_, determines its asymptotic decay and 
defines the power-law exponent of the process. The crit- 
ical value of y that separates these two regimes is 



Vc 



Q+ - 7 

7 — q;_ 



Therefore, when the jumps are exponentially distributed, 
we find a power- law decay for py (y) for any relative value 
of the parameters. 

To test the goodness of the previous analytical deriva- 
tions we have resorted to numerical methods. We have 
simulated the evolution of 1,000,000 different realiza- 
tions of the process ^(t): the relevant PDF is then de- 
rived based on the relative frequency of the outcomes. 
We plot this PDF for r = 100, a quantity that is much 
larger than any characteristic timcscalc of the process, 
since we have set F = 2.0, A = 1.0, A = 2.0, and 7 = 1.0, 
in arbitrary units. For such large values of time, the pro- 
cess is to all effects within the stationary regime, as Fig. [2] 
confirms. 

Figure[2]also shows the coexistence of the two different 
exponents in the decay rate of the stationary PDF, a+ 
and a-. With the choice of the parameters reported 
above a+ = 2.0, a_ = 0.5, and j/c = V^- 



B. Some remarkable general properties 

We have seen in the previous section that if the jump 
density is exponentially distributed the system tends for 
long times to a stationary state where the stationary den- 
sity has a power law decay. We now show that this be- 
havior is not restricted to this particular case but occurs 
under fairly general conditions and election of jump den- 
sity hi'). 

To this end note that in the Laplace domain p(w, t; xq) 
is a well behaved function of lo provided /i(a;) presents 
no pathologies. Nevertheless, /i(w) — >■ as w — >■ 00 and 
hence the inversion to physical space must be done in a 
generalized distributional sense that requires some care: 
indeed the last term of Eq. has no inverse Laplace 
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Figure 2. (Color online) Power-law decay inpy(y). The solid 
(red) line corresponds to the stationary propagator, Eq. (|24|l . 
whereas the dashed (green) line and the dotted (blue) line 
correspond to a decay rate driven exclusively by ot+ or q_ 
respectively. Superimposed to the exact result, and in a good 
agreement with it, we can find (circles) a numerical estima- 
tion of the stationary PDF obtained through Monte Carlo 
simulation. 



transform for x > Xi = -\- Vt and a delta contribution 
appears at that point. This anomaly reflects that there is 
a finite probability that the process arrives to xi before 
any reset or jump has occurred: 

P{X(t) = xil X(0) = xo} = e-(^+^)". (26) 

When r = there is, in addition, a finite chance of find- 
ing the process at a; = after time r, irrespective of the 
initial point. This probability has its origin in the first 
term of Eq. , and amounts 



P{X(t)=0|X(0)=xo} = 



A 



A + A 



f - e 



-(A+A)r 



, (27) 

Remember how these two delta terms have previously 
appeared in Eq. (|20p. Note also that when xq = F = 
the probabilities in Eqs. (|26|) and ((27)) will simply add 
together thus yielding 



>{X(r)=0|X(0) = 0} = 



A + A ■ 



It is also remarkable that the introduction of resets 
makes the system recurrent and ergodic with a stationary 
density that follows from equations and ([T51) : 



p{uj) = lim p(w,t;xo) 
A 



lim sp(cj, s; xq) 

s— ^0 



A 



1 - h{uj) 



(28) 



This ergodicity could be expected on physical grounds 
since the incorporation of this reset mechanism guaran- 
tees that the system will not be driven too far off from the 
origin, no matter how strong the drift and jump mecha- 
nisms are. We elaborate on these properties in App. |B] 



Note that while the limit does not depend on the initial 
condition xq, as expected, it still depends on the drift ve- 
locity F, a fact that is not obvious from Eq. ([TO)) . where 
xq and F always appear side by side. 

In particular, we can use expression (j28p to compute 
the moments of the stationary distribution for a general 
choice of h{-). Thus, if we call the nth-order moment 
of the jump-size PDF h{-) we will have, for instance. 



lim E[A:(r)] 



A/ii 



A 



r Fryr^2l A 2(F + AAti)' 
Inn E A (r) = —^2 + 

r— foo A 



A2 

where E[-] denotes the expectation of its argument. 

We are particularly interested in studying the asymp- 
totic behavior oipx{x) for large value of their arguments. 
By virtue of the Tauberian theorem this amounts to 
studying the limit of Eq. ([^5)) when w — >■ and ft,(w) — > 1. 
Depending on the magnitude of A different scenarios ap- 
pear. Under the assumption that the reset rate is much 
larger than the jump rate, A ^ A, and the velocity^ 
A > F, we expand Eq. (gHl) as 



I 



p(u) 



1 - h(uj) 



. AT A~. . 
1 - - - + -/i(w) 



A A 



A 



This expression leads to 

p{x) 



A 



h{x), 



where we have rejected the terms that do not contribute 
to the asymptotic behavior. 

More compelling is the behavior in the opposite case 
when the above relations are inverted. Concretely, if the 
jump density has finite mean /ii then for small values of 
oj we can approximate 

h{ijj) ^ 1 — niuj + o{ui), 

and hence Eq. ([25)) reads: 

A 



[A^i -I- FJ w + A 
Assume in addition that the ratio 

^ 0(0;), 



A^ti + F 



(29) 



(30) 



i.e., it is at least as small as w is. In this case the Laplace 
inversion yields for large x 

Pi^) - J^^-^'- (31) 

It is interesting to observe the emergence of different 
regimes in pviu) as well, as we increase the reset param- 
eter A. Setting A = (i.e., if we freeze the reset mech- 
anism) the stationary density disappears as befits an in- 
creasing unbounded process. When we introduce a small 
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reset parameter A, the stationary density pviu) exists 
and if condition (|30p is satisfied it will show a heavy-tail 
decay 



Pviy) 



with power- law exponent (3, 



,1+13 ■ 



/3 



A 



A/ii + r 



(32) 



(33) 



Note that this proves that under the above assumptions 
the stochastic process Y{t) is ergodic with a stationary 
density that shows a power-law decay. The critical expo- 
nent depends on a sophisticated combination of all the 
parameters that characterizes the model: the jump and 
reset rates, the drift velocity and the typical jump sizes. 
The result is similar to that reported in Q ■ 

As one keeps increasing the reset rate so that A ^ A 
and A ^ r, the stationary density is dominated by the 
functional form of the jump density, which may or may 
not decay rationally: 

Pviy) ^ -h{\ogy) . 
V 

Finally, for ever larger values of A the pure density part 
is dimmed while the delta term takes over and dominates 
the stationary probability distribution, a result which is 
easy to understand from a physical point of view. 

We finish this discussion by noting that, in the given 
ranges, these expressions reduce to the results obtained 
in the previous section: one just has to recall that from 
Eq. one has fii =7^^. 



IV. EXTREME-EVENT STATISTICS 

In this Section we analyze some statistical properties 
of extreme events of the process X{t). 

We are concerned with the survival probability (SP) of 
the process, 'P[a,b]{T]XQ), 



%,b]{T-xo) = P{a < X{f) < 6,f < t|X(0) - 



2^0} : 



namely, the probability that the process which is initially 
in does not leave the region [a, b]^ a < xq < h, before 
time T. In other words, if we denote by T the first time 
the process exits [a,b] starting from xq, 

T = min {r : AT (r) ^ [a, h] |X(0) = a;o} , (35) 

then the SP is simply the probability that T > t, that is 



'{T <T} = l-r[a,b]{T;Xo) 



(36) 



A related magnitude of interest is the mean exit time 
(MET), the expected value of T : 

T[,,fc](xo) = E[7l. ^ (37) 

Even though these two magnitudes are meaningful for 
any choice of the interval [a, b] we are primarily interested 
in the case a = 0, & > 0: By construction the process 
will tend to reach the upper boundary b steadily, but 
occasionally it returns to the lower boundary a = 0. 

To ease the notation in the sequel we drop the subscript 
in xo, and assume that the process is initially at x. In 
order to get an equation for 7^(0, 6] (t; x) we will resort to 
renewal arguments, very similar to those enumerated in 
Sec. Illli but with the additional requirement that the 
process does not cross the upper boundary at any time. 
With the same notation as in Sec. Illll we will explore the 
three non-overlapping situations : 



(34) 



(i) There is neither a jump nor a reset in the time in- 
terval T, r' > T and r" > r, and the elapsed time is 
not long enough to reach the boundary by the effect 
of the drift, t < {b — x) /T. In this case the process 
will survive with certainty. 

(ii) There is at least one reset event at r' before the first 
jump, t' < t, t' < t" . If t' > {b — x)/r the process 
does not survive, otherwise the survival probability 
is ^[o,6](t-t';0). 

(iii) There is at least one jump before the first reset at 
r", t" < t', t" < t. After this event there is sur- 
vival probability V[ofi]{T — r"; x + Ft' + u) only if 
t" < (b — x) /T and the jump size u is smaller than 
b-x-Tr". 

These three possible contingencies lead to the following 
integral equation: 
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"^[0,6] (t; a;) 



dT"Xe-^^"9{b-x-TT) 



dT"Xe-^^" 



'(A+A)t/ 



J dr'Ae-^^' J dT"\e-^^"V[o,b]{r-T';O)0{b-x-rT') 

duh{u)r[o.b] (r - r"; X + Ft" + u)e{b -x- Tt"), 

e{b-x- Ft) + f dT'Ae''^^+^^^'V[o,b] (r - r'; 0)e{b -x- Ft') 
Jo 

t";x + Ft" + u)e{b-x~FT"). (38) 



+ / dr"Ae-(^+^)"" / duh{u)V[o,b]{r 
Jo Jo 



Next we will consider the Laplace transform of the SP 
with respect to the time variable, 



'P[o.b]is;x) = / dT'P[ofi]{T;x)e' 
Jo 



(39) 



and obtain 



^[0,6] (s; a;) = 



l + AT'io^bjlsiO) 



A + A + s 

f 



1 -e" 



- ^+^+° (b-a:-z) 



X / duh{u)V[o.b] {s; b — z + u) . 



(40) 



Note that in this case we cannot perform the Laplace 
transform with respect to a; in a straightforward way be- 
cause X is restricted to the interval [0, b]. In the following 
section we show how to overcome this limitation. 

We also remind (see App. |D|) that setting s = in 
Eq. (j40p the corresponding integral equation for the mean 
exit time of the process out of the interval, T[o_t](a;), is 
directly obtained: 



Tr[o,6](a;) 



1 + A%bj(0) 
A + A 



1 



fc-2 



dze 



-x — z) 



X / duh{u)T[Qb]ib — z + u). 
Jo 



(41) 



r 



find the general solution of the allied integral equation 



J"(s; z) 



1 + AAjs) 
A + A + s 
A 
f 



1 - e' 



X / duh{u)J-{s;z'~u), 



(42) 



for 2 > 0, with A{s) an arbitrary function of s.0 Notice 
that one recovers "Pio.hj (s; x) via V[o,b] {s;x) = J^(s; b~ x) 
and requiring that A{s) = 7^[o,b] (s; 0). The solution of 
J^(s; z) for values larger than 6, i.e., when a; < 0, is math- 
ematically meaningful but has no physical significance. 

Expression (|42|) is now well suited to be Laplace trans- 
formed in the z variable as well. 



T{s;uj)— / dzT{s;z)e 
Jo 

yielding in this way 

~ 1 1 + Ai(s) 

T{s;uj) = - 



oj \ + A + Fuj + s 
A 



h{iu)T{s;uj). (43) 



A + A + To; + ,s 

Then wc can obtain a closed solution of the problem for 
any functional form of h{-) in the Laplace-Laplace do- 
main: 



1 + AAis) 



1 - h{u}) 



A + Fuj + s 



(44) 



A. The general solution 



Here we consider the solution to Eq. (|40)) with arbitrary 
choice of the jump size PDF h{-). To this end we have to 



B. Exponential jumps 

Let us exemplify the general result in Eq. with the 
exponential case, introduced earlier in Eq. ([T| 



^ One can alternatively obtain this equation by means of the same 
kind of arguments that led to Eq. I l4m . see, e.g., [28l. |41||. 



Note that in Eq. II40II s is a parameter and hence to all effects 
yi(s) plays the role of a constant. 
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Let us begin by writing down the explicit expression 
for J-{s; uj), 

j + uj 1 + AA{s) 



uj Tuj'^ + (A + A + s + w + 7 (A + s) 



1 + AA{s) 1 
A(s) uj 



where 



, , A + A + s + r7 , A(s) 
a±[s) = ± 



2r 



2r 



(45) 



(46) 



A(.s) = y(A + A + s + r7) -4r7(A + s). (47) 

Note that the a± defined above in Eq. ((25|) corresponds 
simply to a± = a±{0). 

The w-Laplace inversion of Eq. P5|) is lengthy but 
straightforward, 



J"(s;z) 



1 + AAjs) 

7 — a_(.s) 



«+(g) - 7 
a+{s) 



and then 

^[0,6] (s;a:^) 



a-{s) 

^ 1 + A^[o.,](s;0) 
A{s) 

^ r «+(g) - 7 

|_ a+{s) 

_^ 7 - a-{s) 
Q;„(s) 



(48) 



-a+(s)(6- 



_ g-a-(s)(b-x)^ 



(49) 



The value of V[Q_t,]{s',0) is obtained after demanding self- 
consistency to expression by letting x = 0. We will 
not pursue further our analysis to detail the explicit form 
of 7-'[o,b](s; 0) since Laplace inversion in the s-variable is 
not possible. 

We will concentrate instead in the s = case, which 
leads to the MET, Tjg^b] (a;), as we have previously stated. 
From Eq. (|49|) we get 

l + AT[o,b](0) 



Tr[o,6](a;) 



F X (a+ — Q!„) 

7 — q;_ 



a+ 



(50) 



The value of T[o^6] (0) must satisfy 



T[o,f,](0) 



l + AT[o,fa](0) 
F X [a^ — «_) 

7 ~ «- 

1 - e 



a+ - 7 



a-b\ 



(51) 



that is. 



T[o,bl(0) 



A 



- a_) 6-"+'' + (c 



- 1 



(52) 



which is always finite since 

A 

a+> - > a-. 

We can find in Fig. [3] some representations of T[o,6] (0) 
where we can observe the opposite effect of changing A 
and F. 




Figure 3. (Color online) Sample representations of T[o,i,](0) 
in terms of A for different choices of P. The solid lines 
were obtained by direct representation of the analytical re- 
sult, whereas the points were computed through numerical 
simulation for: F = 1.0, (red) diamonds; F = 2.5, (green) 
triangles; F = 5.0, (blue) circles; and F — 10.0, (magenta) 
squares. The rest of parameters were b — 1.0, A = 1.0, and 
7 = 1.0. 

Now, we combine Eqs. ([50)) and ([5^ to obtain finally 

1 



(A - FQ;_)e-"+'' + {Ta+ - A)e^"-'' 



«+-7 ^ _^-a+(b-.x) 



a+ 
7 — q;_ 
a- 



l-e' 
1 - e-°'-(''-^' 



(53) 



Figure 2] shows how T[g j,] (x) depends on the starting 
point X, for different intensities of the reset mechanism. 
We can observe how when the reset events are sparse 
T[o ;,](a;) decreases gently, almost linearly, as a direct con- 
sequence of the constant drift term. Reciprocally, when 
there is large reset activity, the escape time is nearly in- 
sensitive to the value of x. This means that in this regime 
the main possibility for the process to leave the interval 
[0, b] comes from the jump mechanism. 

Figures [3] and S] also include numerical estimations of 
the MET for a discrete set of values of the free parameter, 
A and x respectively. Once again we have used a Monte 
Carlo method in which 1,000,000 different realizations 
of the process were averaged to produce the reported re- 
sults. Note that in this case, since the estimation of the 
MET is a problem with a built-in stop mechanism, no 
time cut-off is necessary in the simulation of the process. 
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Figure 4. (Color online) Sample representations of T[o,i,](a;) 
for different values of A. The solid lines were obtained by di- 
rect representation of the analytical result, whereas the points 
were computed through numerical simulation for: A = 0.1, 
(red) diamonds; A = 1.0, (green) triangles, A — 10.0, (blue) 
circles; and finally A = 100.0, (magenta) squares. The rest of 
the parameters were b = 1.0, F = 1.0, A = 1.0, and 7 = 1.0. 



We conclude this section by considering some mathe- 
matical limits on Eq. ((53)) . This will allow us to explore 
more deeply the limiting properties exposed in Figs. [3] 
and m Letting A ^ no reset exists and we recover the 
reported result [3^ \^ : 



X 

+ 



(b-x) 



'-{b-x) 



(A + r7)2 

Letting A — >■ 00 we find 



lim T[o,b] {x) = 

A— >-oo 



alb 



(54) 



(55) 



There exists an alternative, more physical derivation of 
this result by noting that when A — 00 the system is con- 
stantly redirected to 0, and exit can only happen when 
a jump of magnitude greater than b occurs. Hence the 
MET cannot depend on initial position x nor in the drift 
parameter. Furthermore one must also have 



lim T[o,fa](a;) = E [t„ - to\Ji, J„_i < b, Jn > b] 

n=l 
00 

= ^nE[r„]P{Ji,..., J„_i <b,Jn>b} 



g7h 

T' 



and the former result is recovered. This constant value 
is depicted in both Figs. [3] and 2] by a horizontal dotted 
(black) line. 

We can also explore the limits involving the magnitude 



of the drift: 



hmT[o,b](a:)=e^''^{+ ^ 



A A 



1 - e 



{b-x) 



and 



lim T[o,b] {x) = 0. 



(56) 
(57) 



V. CONCLUSIONS 

In this paper we have analyzed in detail the main sta- 
tistical magnitudes of a stochastic process that experi- 
ences random reset events whereupon the system returns 
to the initial point. As the basic stochastic process we 
take a monotonous continuous-time random walk with a 
constant drift: the process increases between the reset 
events, either by the effect of the random jumps, or by 
the action of the deterministic drift. We find that the 
system that results from the balance of these opposite 
effects has interesting properties, like the existence of a 
stationary transition probability density function for any 
choice of jump density and drift strength, or the faculty 
of the model to reproduce powcr-law-like behavior for a 
certain range in the parameter space. If the jump density 
is exponential the transition density can be computed ex- 
plicitly. The long time limit yields a heavy-tailed density 
with a simple power-law exponent. Similarly, exit times 
of the model can also be determined. Analytic results are 
found to be in a remarkable agreement with the numeri- 
cal observations obtained using Monte Carlo analysis to 
average over a large number of realizations. 

Finally note that we have concentrated our efforts in 
the case in which both the waiting times between succes- 
sive reset events and consecutive jumps have a Poissonian 
origin, in the same line of the previous work by Manrubia 
and Zanette Q • The inclusion of alternative innovation 
mechanisms is deferred to a future publication. 
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Appendix A: The general propagator 

Here we analyze in more detail how p{x, r; xq) depend 
on xq. This will allow us to solve Eq. ([TT]) and, eventually, 
to recover Eq. ([T4|) . 

Assume first that x < .tq + Ft. In that case the process 
must at least experience one reset before r since other- 
wise it could not have reached x starting from xq. Since 
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the reset event removes any memory on past locations of 
the process, we must have 

p{x,t;xo) = g{x,T), 

independent of xq. When x > xq + Tt there is in addition 
a possibihty that the process goes from xq to x directly, 
with no reset event in the meanwhile, 

p{x, r; xo) = g{x, r) + Q(x, t; xq). 

Therefore, r; xq) corresponds to the possibility that 
the process reaches x starting from xq by the combined 
action of the random jumps and the constant drift. The 
latter contribution will affect the stochastic process in 
a deterministic way, by decreasing the effective distance 
that the process must travel by means of the random 
jumps alone. Since the jumping mechanism is space- 
homogeneous we must have that, for a given time in- 
terval r, Q{x^t;Xo) must depend on x and xq only via 
the precise combination of variables x — xq — Ft: 

Q{x, r; Xo) = q{x - Xq ~ Tt, t). 

Thus, summing up, p(a;,T;xo) must have the structure 

p{x, r; Xo) = g{x, r) + q{x - .To - Tt, t)0{x - xq - Tt), 

(Al) 

for some functions g{z, r) and q{z, r). Note that the full 
process with drift cannot be recovered from the driftless 
process by the Galilean transformation x — >■ x — Tt. 

We can still gain more insight on g{x, r). This function 
although related with p(x, r; 0) is not the same object. In 
g(x, t) at least one reset event has occurred within the 
time interval r. This means that we know that the time 
of the last reset event is positive, < T*(r) < t, and 
from that point the process has to reach x. Therefore, 
if we define r* = r — T*(t), < r* < r, it follows 
with a similar reasoning to that which led to Eq. (|A1[) 
that we must demand g{x, r) to conform to the following 
structure: 



which reads 



f dT*K{x-TT*,T*)9{x-TT*). 
Jo 



Here the kernel /C(x — Ft*, t*) includes the probability of 
reaching x from the origin after time t* as well as the fact 
that the last reset event was at T* (r) = t — t*; the step 
function simply ensures that x is still accessible. Based 
on the Poissonian nature of the reset events we could 
argue that indeed K{x — Tt*,t*) = Aq{x ~ Ft*, r*), but 
that is not a necessary ingredient at this point: we will 
recover this result as a byproduct of the problem-solving 
procedure. 

This procedure begins with the computation of the 
double Laplace transform of our ansatz for p(x, t; xo), 

p(x,t;xo)= / dT*K{x-rT*,T*)d{x-rT*) 
Jo 

+ q{x - Xo - Ft, r)6l(x - xo - Ft), (A2) 



s; Xo) = -K(uj, s + Flu) 
s 

+ q{uj, s + Fw)e^' 



(A3) 



Inserting this result into the right-hand side of Eq. (jlip 
we obtain 



p(w,s;xo) = 



1 



A + A + s + Fcj 



A 



A + A + s 
A 



1 ~ 



—K{io, s + Tuj) + q{ui, s + Vui) 



A + A + s + Fcj 



h{io)q{uj, s + Fa;)e 



A 



A + A + s s 



-K{uj,s + ruj), 



(A4) 



where 



h{uj) = / duh{u)e-'^'' . 
Jo 

Let us compare the output in Eq. (|A4p with Eq. (|A3P : 
On the one hand equating the terms that contains the 
factor e~"^" we obtain 



q{uj, s + Tuj) 



1 



A + A + .s + Fw 
A 

A + A + s + Fw 



h{Lj)q{uj^ s + Fa;), 



thus we find 



g(w,s) 



1 



A 



1 - h{io) 



+ A + S 



(A5) 



(A6) 



On the other hand considering the terms which are inde- 
pendent from Xo we get 



1 = . ^ 1 
-A' w, s + Fw = 

s ^ ' A+A+s 



A + A 



A'(w, s + Fw) 



and hence 



+ Aq{uj, s + Tuj) 



K{uj, s) ~ Aq{uj, s), 



(A7) 



(A8) 



as we have already anticipated. Therefore, in the Laplace 
domain the solution of the problem reads: 



p(cj,s;xo) 



A 

- +e" 
s 



1 



J A 



1 - h{Lu) 



A + s + Tuj 



Let us analyze the structure of the above equation. 
Note how it can be expressed as: 

p(a;,s;xo) = — po(w,s + A;0) + po(w, s -f A; xo), (A9) 
s 
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with 



Po(w, s; xo) = lim p{u;, s; xq) = — 
A^o ^ 



1 - h{u) 



Tuj- 



Thc interpretation oipo{x, r; xq) is direct: it is the transi- 
tion probabihty of the process between xq and x when the 
reset mechanism is disconnected. Therefore, by appeal- 
ing to well-known properties of Laplace transforms — see 
Eqs. ((C2|) and ([ct)) below—, expression (Kd\ leads to 
Eq. (dH): 

p{x,t;xo)^ / dT*Ae^^^ po{x,T*]0) 

+ PQ{x,T]Xa)e^^'^. 

Note how we have used r* as the integration variable. 
This fact is not incidental at all: r* was defined above 
as the time period comprised between the last reset and 
T, and due to the Poissonian nature of the reset events, 
this interval is exponentially distributed with intensity 
A. The interpretation of Eq. (jl4p is the following: the 
probability of going from Xq to x is the sum of that prob- 
ability when no reset has taken place, plus the probability 
of reaching x after the last reset event. This probabilis- 
tic construction is valid because the propagator does not 
depend on the whole path of the process. 



Appendix B: Markovianity, invariance and 
irreducibility 

This appendix contains the proofs of several probabilis- 
tic properties which have been used in the text. We begin 
with the Markov property. It is well known [l|, Q that a 
time-homogeneous stochastic process X{t) is a Markov 
process if and only if its transition probability satisfies 
the Chapman-Kolmogorov equation (CKE): 



p{x,T + t;xo) 



dxp{x,T; x)p{x,t; Xq), (B1) 



for any f, r. We also note that if X{t) is Markov then 
so it is Y{t). 

A simple way to prove (|B1[) is to use the fact that 
the reset-free process is Markovian and that propaga- 
tors of the process with and without the reset mecha- 
nism, namely p{x,T;xo) and po{x,T; xq), are related by 
Eq. (ITil) . Indeed it implies at once that the right hand 
side of Eq. (jBl[) can be written as a sum of several factors: 



Here 



p{x, T + f; Xo) = Fi + F2 + F3 + F4. 



F2 = / dxpoix,T]x)e-'^^ dT*Ae-'^'' poix,T*]0) 
Jo Jo 



dr*Ae^'^(^+^*)po(a;,T + T*;0) 



dT*Ae-^'' Po{x,T*;0), 



(B2) 



nOO pT 

F3EE dx dT*Ae'^'^po{x,T*;0)po{x,f;xo)e~ 
Jo Jo 



dT*Ae-^(^+^*)po(a;,T*;0), 



(B3) 



and 



poo pr 

Fi= I dx I dT*Ae~^^'pQ{x,T*\Q) 
Jo Jo 

X I df*Ae-^^'po{x,f*;0) 
Jo 

= f dT*Ae-^^'po{x,T*;0) f df*Ae- 

Jo Jo 

dT*Ae-^^'poix,T*;0) (l - e^^H . (B4) 



-At* 





Hence these three factors add to 



F2 + F3 + Fi = J^ dT*Ae-^^'po{x,T*;0) 



+ I dT*Ae-^^'pQ{x,T*;Q) 



dT*Ae-'^^ Po{x,T*;0). (B5) 
We now consider the remaining term 

/•OO 

Fi= dxpo{x,f;xo)po{x,f\xo)e^^^^^''^ 



■ Po{x,T + T;xo)e' 



-A{t+t) 



(B6) 



where we have used that the process without resets is 
Markovian and hence po{x,T; xq) must satisfy its own 
CKE similar to (IbD): 



Po{x,T + t;xo) ^ / dxpoix,T;x)poix,T;xo). (B7) 
Jo 

Finally, using again Eq. (fT?|) we have 

Fi+F2 + F3 + Fi= poix, T + f; xo)e-^^^+^'> 

+ f dT*Ae-^^'po{x,T*;0) 
Jo 

^ p{x,T + f;xo)- (B8) 



Equation (|Bip now follows. 

Next, by letting f ^ 00, Eq. (jBip implies that 



p{x) — / dxp{x,T;x)p{x), 
Jo 

where we remind that p{x) = lYmr^ooPix,T;xo), is in- 
dependent of xo, the initial state. This means that if we 
start with an initial distribution p(-) then for all time 
the density of the process is also p{-). Thus the limit 
distribution is invariant under the dynamics. 

Finally, we prove that the system is irreducible, 
namely, that there is positive probability to visit any re- 
gion A = [x, X + e], X > 0, starting from xq > 0. An 
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analysis of the sample path trajectories along with the 
fact that Jn has a continuous density function h{-) will 
convince the reader that this is indeed the case. For ex- 
ample, \i xq < X then X[t) can reach A with just one 
jump. Consequently the probability P{A,t]Xq) to visit 
A before time t satisfies the lower bound 

px+t—xo 



P{A,T-xo)>e-''Ul 



,-At 



h{x)dx > 0. 



The case x + e < Xq can be handled similarly with the 
only proviso that a reset must occur before r whereupon 
A may be visited via jumps: 

A 



P{A,t)> 



(1- 



-,-At 



A 



A- A 

x-\-e 

h{x)dx > 0. 



A- A 



'1-e 



-At 



Appendix C: Derivation of Eq. ((20)) 

This appendix shows the how the basic result Eq. ([^0)1 
can be recovered by Laplace inversion of Eq. ([T^ . 

We first introduce some notational conventions and list 
two basic properties of the Laplace transform that arc 
used below. For the sake of brevity, we denote by £ and 
the direct and inverse Laplace transform operator 
respectively: 

/(,s) = £[/(T),r,,s], 
/(r) = £-![/>), s,t]. 
The first property is the so-call frequency shifting 

£[/(T)e-^°^T,s] =/> + so), (CI) 
which can be alternatively expressed as 



c- 



f{s + So),S,T 



(C2) 



The second property concerns the interplay between in- 
tegration and Laplace transform: 



f{T)dT,T,S 



or 



£-1 



/(s),s,r 



(C3) 
(C4) 



Thus, we want to compute 

p{x,t;xo) EE C~'^[p{x,s;xo),s,t] 



for 



p{x,s;xo) 



1 



A + A + s 
A 7A 



— 6(x) + S(x — xq) 
s 

(A+3)7X 



s {X + A + sY 
7A 

(A + A + .s)' 



g x+A+s 



Since the Laplace operator is linear, we can write with 
obvious notation 

p{x, r; Xq) ~ Li + L2 + L-i + L4. (C5) 

Let us begin with the simplest term: 

L2 ^ £-1 ^ 



A + A + s 



s,T 6{x-xo) ^e~'^^+^'>^d{x~XQ), 

(C6) 

where we have used the property (jC2p for sq = A + A and 
the fact that 



£-1 



1 

-,S,T 



1. 



(C7) 



The next term Li is very similar to L2 using partial 
fraction decomposition: 

-,S,T 5{x) 



Li = £-1 
A 



A + A 
A 



s A + A + s 
£- ' 



s A + A + s 



, S, T 



Six) 



1 - e 



-(A+A)r 



A + A 
The fourth term reads 



Six) 



(C8) 



£-1 



7A 



(A + A + s)' 



_g- A+A+. T(a:-3:o)^ 



[X - Xo) 



= 7Ae'(^+A)^£- 
x£-^l 



9ix - Xo) 



(x - Xo) 



,-(A+A)t-7(x-xo) 



7Ar 

X — Xq 



X /i (2v/Ar7(x - xo)) 0(x - xq). (C9) 
where we have used again Eq. ()C2p and also that: 

1 c 



£-1 



re=,s,r 



/i (2VH7) , (CIO) 



see [43. Here /i(-) is the modified Besscl function of 
order 1. 

Only the term remains. It can be handled by using 
Eqs. dell), (IC2]) and ([CTO| consecutively: 



L. = £-1 



A 7A 

s (A + A + s) 



26 ^+A+s ^ T 



= A7A /" df£-i 
"'0 

= A7A / dfC~^ 



(A + A + s) 

1 A7a; _ 

= ,s,T 



g A+A+s c 7- 



-(A+A)t-72: 



A / ■ dfe-(^+^)^-^y ^/i (2v/A?^) 



(Cll) 



Putting together all the terms, Eq. (PU|) follows. 
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Appendix D: From survival probabilities to mean 
exit times 



According to the results of Sec. iTVl and Eq. ((36l) . the 
mean exit time f,](xo) may be computed as follows 



1*00 1*00 

la,b]i^0)= TdP{T<T} = - Tdr[,,b]{T;Xo) 

Jo Jo 



where in the last equality we have assumed that 
'Pla.b](j"yXo) decays appropriately as r — > 00 to perform 



safely integration by parts. Now, note that by definition 



^[Q,b](s = 0;a;o) = / 'P[a.b]{T; xo)e ^^dr 
Jo 

= / 'P[a.b](^'^^o)dT, 

Jo 

and thus finally 

T[a,b](a;o) = ^[a,6](s = 0;a;o). 



s=0 



(Dl) 
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